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Abstract. The calculation of the trajectories in the Sun-Earth-Mars system will be 
performed in two different models, both in the framework of Newtonian mechanics. 
First model is well-known Copcrnican system, which assumes the Sun is at rest and all 
the planets orbit around it. Second one is less-known model developed by Tycho Brahe 
(1546-1601), according to which the Earth stands still, the Sun orbits around the Earth, 
and other planets orbit around the Sun. The term "Neo-tychonian system" refers to 
the assumption that orbits of distant masses around the Earth are synchronized with 
the Sun's orbit. It is the aim of this paper to show the kincmatical and dynamical 
equivalence of these systems, under the assumption of Mach's principle. 
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1. Introduction 

The discussion of motion of celestial bodies is one of the most interesting episodes in 
the history of science. There are two diametrically opposite schools of thought: one 
that assumes that the Sun stands still, and Earth and other planets orbit around it; and 
another that assumes that the Earth stands still, and Sun and other planets in some 
manner orbit around the Earth. The first school of thought comes from Aristarchus (310- 
230 bc) and is generally addressed as heliocentrism, another from Ptolemy (90-168 bc) 
and is generally known as geocentrism. Since Aristotle, the ultimate authority in science 
for more than two millennia, accepted the geocentric assumption, it became dominant 
viewpoint among scientists of the time. The turnover came with Copernicus (so-called 
"Copernican revolution" ) who in his work De Revolutionibus proposed a hypothesis that 
the Sun stands in the middle of the known Universe, and that Earth orbits around it, 
together with other planets. Copernicus' system was merely better than Ptolemy's, 
because Copernicus assumed the trajectories of the planets are perfect circles, and 
required the same number of epicycles (sometimes even more) as Ptolemy's model [I]. 
The accuracy of Ptolomy's model is still a subject of vivid debates among historians of 
science [2]. 

The next episode in this controversy is Kepler's system with elliptical orbits of 
planets around the Sun. That system did not require epicycles, it was precise and 
elegant. It is therefore general view that Kepler's work finally settled the question 
whether it is the Sun or the Earth that moves. But what is less known is that Tycho 
Brahe, Kepler's tutor, developed a geostatic system that was just as accurate and elegant 
as Kepler's: the Sun orbits around the Earth, and all the other planets orbit around 
the Sun. The trajectories are ellipses, and all the Kepler's laws are satisfied. In that 
moment of history, the Kepler's and Brahe's models were completely equivalent and 
equally elegant, since neither of them could explain the mechanism and reason why the 
orbits are the way they are. It had to wait for Newton. 

Sir Isaac Newton, as it is generally considered, gave ultimate explanation of 
planetary motions that was in accord with Kepler's model, and excluded Brahe's one. 
The laws of motions and the inverse square law of gravity could reproduce all the 
observed data only with the assumption that the Sun (i.e. the center of mass of the 
system, which can be very well approximated by the center of the Sun) stands still, and 
all planets move around it. According to Newton's laws, it is impossible for small Earth 
to keep the big Sun in its orbit: the gravitational pull is just too weak. This argument 
is very strong, and it seemed to settle the question for good. 

But in the end of 19th century, the famous physicist and philosopher Ernst Mach 
(1839-1916) came with the principle which states the equivalence of non-inertial frames. 
Using the famous "Newton's bucket" argument, Mach argues that all so-called pseudo- 
forces (forces which result from accelerated motion of the reference frame) are in fact 
real forces originating form the accelerated motion of distant masses in the Universe, as 
observed by the observer in the non-inertial frame. Some go even further, stating that 
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"every single physical property and behavioral aspect of isolated systems is determined 
by the whole Universe" [3]. According to Mach's principle, the Earth could be considered 
as the "pivot point" of the Universe: the fact that the Universe is orbiting around the 
Earth will create the exact same forces that we usually ascribe to the motion of the 
Earth. 

Mach's principle played a major role in the development of the Einstein's General 
Theory of Relavity [3], as well as other developments in gravitation theory, and has 
inspired some interesting experiments [5]. This principle still serves as a guideline for 
some physicists who attempt to reformulate ( "Machianize" ) Newtonian dynamics [6], [7] , 
or try to construct new theories of mechanics [8]. Some arguments and critiques against 
Mach's principle have also been raised [9]. Since the time of it's original appearance 
[TOl [TTl [12] , Mach's principle has been reformulated in numbers of different ways [13l [14] . 
For the purpose of this paper, we will only focus on the one of the consequences of Mach's 
principle: that the inertial forces can be seen as resulting from real interactions with 
distant matter in the Universe, as was for example shown by A. Zylbersztajn |15j . 

The only question remains: are these forces by themselves enough to explain all 
translational motions that we observe from Earth, and can they reproduce the Tycho 
Brahe's model? The discussion in this paper will show that the answer to this question 
is positive. In order to demonstrate it, we will consider the Sun-Earth-Mars system. 

The paper is organized as follows. In section [2] an overview of two-body problem in 
the central potential and of Kepler's problem is given. In section [3] the calculations 
of Earth's and Mars' trajectories are performed in the heliocentric system, both 
analytically (by applying the results from previous section) and numerically. In section 
[4] the calculations of Sun's and Mars' trajectories are performed in geocentric system, 
due to the presence of pseudo-potential originating from the fact of accelerated motion 
of the Universe. Finally, the conclusion of the analysis is given. 



2. Two-body problem in the central potential 

2.1. General overview 

We start with the overview of two body problem in Newtonian mechanics. Although 
there are alternative and simpler ways to solve this problem (T6J [17], we will follow the 
usual textbook approach [18j [19]. The Lagrangian of the system reads 

L = l - mi il + X -m 2 i\ - t/(|r x - r 2 |) , (2.1) 

where U is potential energy that depends only on the magnitude of the difference of 
radii vectors (so-called central potential). We can easily rewrite this equation in terms 
of relative position vector r = ri — r 2 , and let the origin be at the centre of mass, i.e. 
mil - ! + m 2 r 2 = 0. Solution of these equations are 

m 2 mi 
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The Lagrangian (12. ip so becomes 



L = - U(r) , (2.3) 

where r = |r| and \x is the reduced mass, 

111 , , 

- = — + — 2.4 

H mi m 2 

In that manner, the two-body problem is reduced to one-body problem of particle with 
coordinate r and mass [i in the potential U(r). 

Using polar coordinates, the Lagrangian f]2 .31) can be written as: 

L = l/i(r 2 + r 2 2 ) - U(r) (2.5) 

One can immediately notice that variable is cyclic (it does not appear in the 
Lagrangian explicitly). Consequence of that fact is momentum conservation law, since 
(d/dt)(dL/d(f)) = dL/d<j> = 0. Therefore, 
dL 

£ = — - = ur 2 d> = const. (2.6) 

d<f) K J 

is the integral of motion. 

In order to find a solution for the trajectory of a particle, it is not necessary to 
explicitly write down the Euler-Lagrange equations. Instead, one can use the energy 
conservation law, 

E = i/i(r 2 + r 2 4> 2 ) + U{r) = ^f 2 + ^ + U(r) (2.7) 

Straightforward integration of (12 .7p gives equation for the trajectory, 

£ dr/r 2 



[E - U(r)] - £ 2 /r 2 



{2.1 



2.2. Kepler's problem 



Let us now consider the particle in the potential 

U(r) = --, (2.9) 
r 

generally known as Kepler's problem. Since our primary interest is in the planetary 
motions under the influence of gravity, we will take k > 0. Integration of eq. (I2.8P for 
that potential gives: 

- = l + ecos0, (2.10) 
r 

where 2p is called lactus rectum of the orbit, and e is eccentricity. These quantities are 
given by 



/■2 



2E£ 2 

* = W e = f + w (2 ' n) 
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Expression (I2.10p is the equation of a conic section with one focus in the origin. For 
E < and e < 1 the orbit is an ellipse. 

One can also determine minimal and maximal distances from the source of the 
potential, called perihelion and aphelion, respectively: 

Train Z ' j Tmax Z • (2.12) 

1 + e 1 — e 

These parameters can be directly observed, and often are used to test a model or a 
theory regarding planetary motions. 

3. Earth and Mars in heliocentric perspective 

According to Newton's law of gravity, the force between two massive objects reads: 

Gm x m 2 

F = ~, ^(ri-r 2 ). (3.1) 

|ri - r 2 | 

Which leads to a potential (F = —VU) 

U{ \ tl -r,\) = -^2SSl (3.2) 
|ri - r 2 | 

This is obviously Kepler's potential ( 12. 9 p with k = Gmim2, where G is Newton's 
gravitational constant. 

Since the Sun is more than 5 orders of magnitude more massive than Earth and 
Mars, we will in all future analysis use the approximation 

fittrrii, (3.3) 

where m ; is mass of the observed planet. For the same reason, gravitational interaction 
between Earth and Mars can be neglected, since it is negligible compared with the 
interaction between Earth/Mars and the Sun. 

Using these assumptions, we can write down corresponding Lagrangians, 

1 . 2 Gm E M s 
Les = ^ m E^Es + , 

2 r E s 

t 1 -2 . Gm M M s 

Lms = zz m M^ ms H > ( 3 - 4 ) 

2 r M s 

where m E and tum are masses of Earth and Mars, respectively. Subscripts ES (MS) 
correspond to the motion of Earth (Mars) with respect to the Sun. These trajectories 
can be calculated using the exact solution ( 12.101) with appropriate strength constants 
k and initial conditions which determine E and £. Another way is to solve the Euler- 
Lagrange equations numerically, using astronomical parameters [20] (e.g. aphelion and 
perihelion of Earth/Mars) to choose the inital conditions that fit the observed data. 
The former has been done using Wolfram Mathematica package. The result is shown 
on Fig. [TJ 

For latter comparison, one could write out the expressions for the e and p parameters 
for the Earth. Putting the expressions for energy (12.71) and momentum (12. 6p into 
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Figure 1. Trajectories of Earth and Mars in heliocentric system over the period of 
2 years. Blue and red lines represent Earth's and Mars' orbits, respectively (color 
online) . 



Equations (12. lip it is straightforward to obtain 

2 r 4 



V 



GMc 



2GM S (j) 2 r 3 - r 2 2 r 4 - ^r 6 
Y GHfl ' (3 - 5) 

where 0, f and r are angular velocity, radial velocity and distance respectively, taken in 
the same moment of time (e.g. in t = 0). 

Fig. [2] displays motion of the Mars as viewed from the Earth, gained by trivial 
coordinate transformation 

rjifs(t) = -T E s(t) + r M s{t) , (3-6) 

where r ES (t) and r M s(t) are solutions of Euler-Lagrange equations for the Lagrangians 
(13. 4p . Equation (13. 6p is just the mathematical expression of the Tycho Brahe's claim. 
Retrograde motion of the Mars can be useful in the attempt to understand and determine 
orbital parameters, as was qualitative and quantitative shown by B. Thompson [2T] . 

The acceleration that Earth experiences due to the gravitational force of the Sun 
is usually referred as centripetal acceleration and is given by 
F cp GM S „ 

a c P — — 5 — r ES > (3-7) 

m E r ES 
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Figure 2. Trajectory of the Mars as seen from the Earth over the period of 7 years. 
Calculation of this trajectory is done numerically in the heliocentric system. 

where r is the unit vector in the direction of vector r, r_Es(t) is radius vector describing 
motion of Earth with respect to the Sun, and F cp is centripetal force, i.e. the force that 
causes the motion. 

4. Sun and Mars in geocentric perspective 

4-1. The pseudo-potential 

From the heliocentric perspective, the fact that the Earth moves around the Sun results 
with centrifugal pseudo-force, observed only by the observer on the Earth. But if we 
apply Mach's principle to the geocentric viewpoint, one is obliged to speak about the real 
forces resulting from the fact that the Universe as a whole moves around the observer 
sitting on the stationary Earth. Although these forces will further be considered as the 
real forces, we well keep the usual terminology and call them pseudo-forces, for the sake 
of convenience. Our focus here will be on the annual orbits, not on diurnal rotation 
which requires some additional physical assumptions [8] [22] that are beyond the scope 
of this paper. 

The Universe is regarded as an (N + l)-particle system (N celestial bodies plus 
planet Earth). From the point of a stationary Earth, one can write down the Lagrangian 
that describes the motions of celestial bodies: 

Z i=l Z i=l 'ij i=l 'i 
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where = jr^ — r ^ | , U ps stands for pseudo-potential, satisfying F ps = — VU ps . F ps is 
the pseudo-force given by 

N 

F ps = -m a cp ,i , (4.2) 
i=i 

where a cPji is centripetal acceleration for given celestial body (with respect to the Earth) 
and m is a mass of the object that is subjected to this force. It's easy to notice that the 
dominant contribution in these sums comes from the Sun. The close objects (planets, 
moons, etc) are much less massive than the Sun, and massive object are much further 
away. The same approximation is implicitly used in section |3j 

In the Machian picture, the centripetal acceleration is a mere relative quantity, 
describing the rate of change of relative velocity. Therefore, centripetal acceleration of 
the Sun with respect to Earth is given by Equation (13. 7p . with r ES = —r$E- All that 
considered, Equation (I4.2[) becomes 

Fps = _£^ fsE , (4.3) 

r SE 

where rsE{t) describes the motion of the Sun around the Earth, and m is the mass of 
the body under consideration. 

We can now finally write down the pseudo-potential which influences every body 
observed by still observer on Earth: 

U ps (r) = ^^r SE -r, (4.4) 

r SE 

where r(t) describes motion of particle of mass m with respect to the Earth. Notice 
that this is not a central potential. 



4-2. Sun in Earth's pseudo-potential 

In order to determine Sun's orbit in Earth's pseudo-potential, one needs to take 
dominant contributions of the Lagrangian (14.11) . as was explained earlier. Taking into 
account the expression for pseudo-potential given in Equation ( I4.4p . one ends up with 

1 „ GMl , . 
L SE = -M s r 2 SE s - . 4.5 

2 r SE 

This Lagrangian has the exact same form as the reduced Lagrangian (12.31) . That 
means that we can immediately determine the orbit by means of Equations (12. lip by 
substituting /i = M s and k = GMg. This leads to the following result (subscript SE 
will be omitted): 

2 r 4 
P= GM~s 1 
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which is the exact equivalent as the previous result given in Equations (I3.5p . since 
0, f and r are relative quantities, by definition equivalent in both models. We can 
therefore conclude that the Sun's orbit in the Earth's pseudo-potential is equivalent to 
that observed from the Earth in the heliocentric system. 
It remains to show the same thing for Mars' orbit. 



4-3. Mars in Earth's pseudo-potential 

In the similar way as before, we take the dominant contributions of Lagrangian (14. ip 
together with Equation (14. 4 j) and form the following Lagrangian 



'ME 



-m M r ME 



Gm M M s 



Gm M M s 



-r SE -r ME , (4.7) 

\*me - *se\ r- SE 

where subscript ME refers to the motion of Mars with respect to the Earth, and rs#(i) 
is the solution of Euler-Lagrange equations for the Lagrangian (14. 5p . 




Astronomical units 



Figure 3. Trajectories of the Sun (yellow) and the Mars (red) moving in Earth's 
pseudo-potential over the period of 7 years (color online). Calculation of this trajectory 
is performed numerically in the geocentric system. 



The Euler-Lagrange equations for vme($) using Lagrangian ( 14. 7 j) are too 
complicated to be solved analytically, but they can easily be solved numerically. The 
numerical solutions for the equations of motion for both the Sun and Mars are displayed 
in Fig. |3j The equivalence of trajectories gained in two different ways is obvious, 
justifying the model proposed by Tycho Brahe. 
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5. Conclusion 



The analysis of planetary motions has been performed in the Newtonian framework with 
the assumption of Mach's principle. The kinematical equivalence of the Copernican 
(heliocentric) and the Neo-tychonian (geocentric) systems is shown to be a consequence 
of the presence of pseudo-potential ( 14.4p in the geocentric system, which, according 
to Mach, must be regarded as the real potential originating from the fact of the 
simultaneous acceleration of the Universe. This analysis can be done on any other 
celestial body observed from the Earth. Since Sun and Mars are chosen arbitrarily, and 
there is nothing special about Mars, one can expect to come up with the same general 
conclusion. 

There is another interesting remark that follows from this analysis. If one could 
put the whole Universe in accelerated motion around the Earth, the pseudo-potential 
corresponding to pseudo-force (14. 2 p will immediately be generated. That same pseudo- 
potential then causes the Universe to stay in that very state of motion, without any 
need of exterior forces acting on it. 



Acknowledgments 



This work is supported by Ministry of Science, Sports and Technology under contract 
number 119-0982930-1016. 



References 



[1] Koestler A 1959 The Sleepwalkers: A History of Man's Changing Vision of the Universe 

(Hutchinson: MacMillan), pp. 194-195 
[2] Rawlins D 1987 "Ancient hclioccntrists, Ptolemy, and the equant" Am. J. Phys. 55 235-239 
[3] Rosen J. 1981 "Extended Mach principle" Am. J. Phys. 49 258-264 

[4] Newburgh R 2007 "Inertial forces, absolute space, and Mach's principle: The genesis of relativity" 

Am. J. Phys. 75 427-430 
[5] Lichtenegger H and Mashhoon B 2004 "Mach's principle" arXiv:physics/0407078 [physics. hist- ph] 
[6] Hood C G 1970 "A Reformulation of Newtonian Dynamics" Am. J. Phys. 38 438-442 
[7] Barbour J 1974 "Relative-distance Machian theories" Nature 249 328-329 
[8] Assis A K T 1999 Relational Mechanics (Montreal: Aperion) 

[9] Hartman H I and Nissim-Sabat C 2003 "On Mach's critique of Newton and Copernicus" Am. J. 
Phys. 71 1163-1168 

[10] Mach E 1872 Die Geschichte und die Wurzel des Satzes von der Erhaltung der Arbeit, (Prague: 
Calve) 

[11] Mach E 1883 Die Mechanik in ihrer Entwickelung historisch-kritisch dargestellt, (Leipzig: 
Brockhaus) 

[12] Mach E 1911 History and Root of the Principle of the Conservation of Energy (Chicago: Open 
Court) 

[13] Rovelli C 2004 Quantum Gravity (Cambridge University Press) p 75 

[14] Barbour J 2010 "The definition of Mach's principle" larXiv: 1007.33681 [gr-qc] 

[15] Zylbersztajn A 1994 "Newton's absolute space, Mach's principle and the possible reality of fictitious 

forces", Eur. J. Phys. 15, 1-8 
[16] Hauser W 1985 "On planetary motion" Am. J. Phys. 53 905-907 



Newton- Machian analysis of Neo-tychonian model of planetary motions 



11 



[17] Gauthier N 1986 "Planetary orbits" Am. J. Phys. 54 203 

[18] Landau L D and Lifshiz E M 1976 Mechanics 3rd ed (Oxford: Butterworth-Heinemann) pp. 25-40 
[19] Goldstein H 1980 Classical Mechanics 2nd cd (Addison- Wesley) pp. 70-102 

[20] Handbook of Chemistry and Physics (49th ed) ed Weast R C 1968 (Cleveland Ohio: The Chemical 

Rubber Company), pp. F 145-146 
[21] Thompson B G 2005 "Using retrograde motion to understand and determine orbital parameters" 

Am. J. Phys. 73 1023-1029 
[22] Veto B 2011, "Gravitomagnctic field of the universe and Coriolis force on the rotating Earth" Eur. 

J. Phys. 32 1323-1329 



